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Abstract 

We study the problem of solving integration-by-parts recurrence relations for a 
given class of Feynman integrals which is characterized by an arbitrary polynomial 
in the numerator and arbitrary integer powers of propagators, i.e., the problem of 
expressing any Feynman integral from this class as a linear combination of mas- 
ter integrals. We show how the parametric representation invented by Baikov ^ 
can be used to characterize the master integrals and to construct an algorithm for 
evaluating the corresponding coefficient functions. To illustrate this procedure we 
use simple one-loop examples as well as the class of diagrams appearing in the 
calculation of the two-loop heavy quark potential. 

PACS numbers: 02.70.-c, 12.38.-t, 14.65.-q 



1 Introduction 

In the recent years the art of evaluating multi-loop Feynman integrals has been driven to 
an impressive high level and very sophisticated methods and tools have been developed in 
order to cope with the enormous complexity one encounters at higher orders (see, e.g., the 
recent reviews |2]). A standard practical problem is the evaluation of a class of Feynman 
integrals corresponding to a given graph. The integrals differ by the integer powers of the 
propagators and the polynomials in numerators. Most methods heavily rely on the use 
of recurrence relations derived with the help of integration-by-parts (IBP) |3j identities. 
The recurrence relations are used to express a complicated integral in terms of simpler 
ones and, after repeated use, finally in terms of a small set of integrals, so-called master 
integrals, which can not further be reduced. Only for the latter a complicated integration 
is possibly necessary. 

Any dimensionally regularized Feynman integral corresponding to a given graph, with 
the space-time dimension d = 4 — 2e as a regularization parameter , can be represented 
in the form 

^ ^ f d% . . . d'h 

where ki, i = 1, . . . , h, are loop momenta, rii are integer indices, underlined letters denote 
multi-indices, i.e., n = (ni, . . . , tin), and the denominators are given by 

Da= Y: A^^P^■p,-ml (2) 

with a = 1, . . . ,N. Irreducible numerators which cannot be linearly expressed through 
the given set of the denominators Da of the propagators are naturally treated as extra 
denominators raised to a negative power. The momenta pi are either the loop momenta 
Pi = ki, i = 1, . . . , h, or external momenta Ph+i, ... of the graph so that the denominators 
can be either quadratic or linear^ in the loop momenta ki. 

The IBP identities are obtained by applying the operator (d/dpi) ■ pj {i = 1, . . . ,h 
and j = 1, . . . , h, . . .) to the integrand in Eq. ([T]). Due to the properties of dimensional 
regularization the resulting integral is zero which can symbolically be written as 

P,,(/+,/-)F(n) = 0. (3) 

Pij are polynomials in the operators that increase and lower indices, and I~, where 

I+F(...,n„...)=F(...,n, + l,...), 

I-F(...,n„,...)=F(...,n,-l,...). (4) 

The basic idea is to use the relations of Eq. Q in order to express any integral with 
a given value of the multi-index n in terms of simpler integrals corresponding to some 

^Linear denominators usually appear in asymptotic expansions of Minkowski space Feynman integrals 
within the strategy of expansion by regions [S]- 
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finite family of multi-indices nj, where, usually, = {rin, . . . , n^Tv) consists^ of riia = 0, or 
1, or a negative integer value. Having in mind the experience collected in the process of 
solving recurrence relations for various classes of Feynman integrals, one may hope that, 
for any given class of Feynman integrals, the recursion problem can be solved, i.e., any 
Feynman integral F{n), with given values of its indices, can be represented as a finite 
linear combination 

F{n) = J2c,{n)F{n^, (5) 

i>l 

with the normalization condition 

CiiRj) = 6ij . (6) 

The integrals F{nj) are called master, basic or irreducible integrals. Our experience tells 
us that the coefficient functions Ci{n) turn out to be rational functions of the dimension d, 
the masses and the kinematical invariants build from the external momenta, so that the 
non-trivial analytical dependence of the Feynman integrals on the dimensionful quantities 
and d is completely contained in the master integrals. 

To solve the problem of the reduction it is necessary to identify the master integrals 
and to construct an algorithm which allows to calculate the coefficients Q(n). In concrete 
situations, the realization of the reduction procedure to a set of master integrals turns 
out to be far from straightforward. Examples of recent attempts to construct systematic 
procedures for solving IBP recurrence relations can be found in [HI HI |H1 E] ■ 

An alternative approach to solve recurrence relations in a systematic way was devel- 
oped in pp. This approach is based on an appropriate integral representation for the 
coefficient functions Ci{n). In the case of vacuum Feynman integrals, this representation 
takes the form pp 

c.(n)-/.../%^[P(x)]^'^-'^-^^/% (7) 

where the parametric integrals have to fulfill the essential condition that IBP can be 
applied. As one can see in the examples discussed in Section El the integrals in Eq. ((7j) 
can be either performed as closed contour integrals in the complex plane, or as iterated 
integrations over real one-dimensional domains. The latter is typically between roots of 
a quadratic polynomial. Furthermore, if the singularity of the Xj-integration is given by 
the factor (n^ = 1,2,...) it is unreasonable to choose this point as a boundary of 
a segment-integration as the corresponding divergence is not regularized. However, such 
a choice turns out to be possible if (probably, after some intermediate integrations) an 
additional factor of the type x^^", where k is integer, appears. 

^Examples are known (see, e.g., Ref. |^) where for master integrals some of the indices are equal to 
two, and not only to zero and one. However, in these situation, one can switch to another set of master 
integrals corresponding to indices equal to zero or a negative integer value. 
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The basic polynomial in Eq. ((Tj) is given by TQ 
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where the determinant is taken with respect to the indices i and j. Thus, P is a homo- 
geneous polynomial of degree h in the variables Xi and the masses. The matrix A'^ in 
Eq. (jSl) is defined as follows: A^^ as introduced in Eq. Q is defined for i > j. Let us then 
consider the quadratic N x N matrix A (with N = h{h + l)/2), where the first index is 
labeled by pairs with i > j, and the second index is a. The corresponding inverse 
matrix {A~^yj (where again i > j) satisfies 



Now A^J is defined for all i and j as the symmetrical extension of [A^^YJ . 

As was demonstrated in jT], the representation ((Tj) satisfies Eq. (j21), provided one can 
use IBP in the parametric integrals. We should remark that the overall normalization 
in Eq. ((7j) is not fixed in advance but is adjusted after the construction of the coefficient 
functions, as a virtue of Eq. (jH)). For the same reason, the basic polynomial is determined 
up to a factor which is independent of the variables Xj. 

General Feynman integrals can be reduced to the vacuum case ^1 QS]- In case an 
external momentum is reduced to a mass shell, = mf (where the mass rrij is one of 
the internal masses of the diagram), one includes into the procedure all the terms of the 
formal Taylor expansion in at pf = mf. In the case of a general external momentum 
squared one can still consider similar terms of the Taylor expansion in pf — Mf where, 
this time. Mi is not equal to some of the internal masses. This Taylor expansion is 
indeed a Taylor series (rather than a formal series) because the point pf = Mf is not 
singular. In this situation, one is usually interested only in the value at pj = Mf and not 
in the derivatives at this general point. Thus for the corresponding Feynman integrals 
one considers only the value rij = 1 of the corresponding index n^. The transition to a 
vacuum problem effectively increases the number of loops which enters the exponent of 
the basic polynomial P in Eq. where the loop number h refers to the corresponding 
vacuum Feynman graph. For example, in a recursive problem for a class of propagator 
diagrams, the transition to the corresponding vacuum problem, which is obtained after 
gluing the external lines, increases the loop number by one. We would like to stress that 
this translation to the level of vacuum Feynman integrals is not necessarily accompanied 
by a corresponding vacuum graph. However, in the language of Feynman integrals defined 
through Eq. ((T)) such a description is always possible. 

It should be mentioned that although Eq. ((Tj) has been applied to solve IBP recurrence 
relations in some physical situations ^] no instructions how to apply this represen- 
tation are available in the literature. 

The goal of our paper is twofold. First we shall describe how the parametric repre- 
sentation of Eq. can be used to characterize the master integrals and to construct an 
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algorithm for evaluating the coefficients of these master integrals. Second, we shall choose 
several typical examples to exemplify the method. As a non-trivial two-loop example, we 
discuss in detail the class of diagrams appearing in the calculation of the two-loop heavy 
quark potential. The main motivation for this is that in Refs. fH^ IT^ no complete pro- 
cedure has been presented but only the necessary Feynman integrals were calculated. A 
reduction procedure is described in the thesis [T3] which is based on IBP accompanied 
by the method developed in jB]. In Ref. where l/(mqr^) corrections to the potential 
were evaluated, a standard approach to solve the IBP relations has been used and an 
algorithm has been developed to evaluate the adequate subclass of the diagrams. In this 
paper we present a complete procedure for evaluating a general diagram of this class which 
is certainly useful if higher dimensional operators are considered within the framework of 
non-relativistic QCD. 

The outline of this paper is as follows: in the next section, we discuss the identification 
of the master integrals and present general prescriptions for applying Eq. ((7j). In Section |2l 
we discuss the one- and two-loop examples and, finally, in Section HI we present our 
conclusions. 

2 Classification of master integrals and construction 
of coefficient functions 

In this section we present practical prescriptions to perform a reduction procedure with 
the help of Eq. ((7j). We shall assume that any Feynman integral can be decomposed 
according to Eqs. (0) and The method decomposes into two parts: the classification 
of the master integrals and the construction of the coefficient functions. 

2.1 Identifying master integrals 

Let us in a ffist step consider the integrals where the indices corresponding to irreducible 
numerators are set to zero. The other indices we allow to be either zero or one. In the 
following we nullify all scaleless integrals (e.g., massless tadpoles) which is natural within 
the framework of dimensional regularization. 

Before analyzing the remaining non-zero integrals with the indices one and zero, let us 
remember that the goal of every known reduction is to reduce some of the indices to zero. 
Indeed, experience tells us that a master integral F{nj) never appears in the reduction 
of a given Feynman integral F{n) if Uj < and Uij > 0. For this reason we adopt the 
following natural condition for the coefficient function Cj(n) of F{nj): 

Ci{n) =0 if rij < and riij = 1 for j = 1, . . . , . (10) 

This condition is easily realized in an automatic way by choosing the integration domain 
for Xj in Eq. ((Tj) as a closed contour around the origin in the complex Xj-plane, which 
from now on will always be implied. As a consequence, the integration over the variables 
Xj in Eq. (jZj), where Uij = 1 in the corresponding master integral, is performed with the 
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help of the Cauchy theorem and thus reduces to Taylor expansions of order rij — 1 of the 
integrand in Xj. Eventually, Eq. (|7j) results in a sum of terms of the form 

f...[mx)Y-^ n (11) 

where z = (c? — /i — l)/2, is a non-negative integer, n'j are some integer exponents, and 
Pi{x) is obtained from P{x) after setting to zero all variables Xj where riij = 1. 

Let us now analyze the remaining candidates in view of being master integrals with 
the help of the representation (fTTj). If Pj = for a given candidate F{nj) {i.e., for 
with riij = 1 or 0) it is naturally to put the corresponding coefficient function to zero. In 
other words, such a Feynman integral is recognized as a reducible integral in the sense 
of the reduction problem. This is the simplest condition of reducibility. A more general 
condition is that any resulting integral of Eq. (fTTj) involves an integral without scale. In 
this situation it is natural to prescribe a zero value for it^ and, therefore, consider the 
given Feynman integral F{n^ as reducible in the sense of the reduction problem and the 
corresponding decomposition given in Eq. It can happen that a scaleless integral 
in Eq. (|7j) can either be seen immediately, i.e., an integral over a "pure" power of some 
variable arises, or only after some intermediate integrations. 

After trying all possible combinations of zero and one we obtain a list of master 
integrals, which fixes the terms on the right-hand side of Eq. (jSj). 

However, this is not all. Sometimes, in addition to a given master integral F{n^ with 
Ui consisting of riij = 1 or rtij = 0, one has to consider master integrals which differ 
from F{n^) by some indices Uij < 0. The number of such additional master integrals is 
dictated by the degree of the polynomial Pi with respect to some of the parameters Xj. 
In practice it turns out to be advantageous to complete the list of the master integrals 
in the second step where algorithms for the computation of the coefficient functions Ci{n) 
are constructed. 

Sometimes it is obvious that a given candidate is reducible due to other arguments like 
the simple application of the triangle rule. However, it turns out to be rather convenient 
to see this property directly by analyzing the corresponding "reduced" polynomial Pi, as 
this provides a useful additional check. 



2.2 Constructing coefficient functions 

The set of the master integrals with the indices zero and one as obtained in the previous 
Subsection is partially ordered in the natural way. We write F{ni) < F{n2) if riij < n2j 
for all j and the strict inequality holds at least for one index. In this way all master 

•^At this point we once again apply a commonly accepted agreement within dimensional regularization 
to put any scaleless integral to zero. Observe that such a prescription, in an extended form, is also 
successfully applied within the strategy of expansion by regions [S] when expanding Feynman integral in 
various limits of masses and momenta. Moreover scaleless integrals that are not dimensionally regularized 
(e.g., connected with the contribution of the potential region) are also consistently put to zero with this 
strategy. 
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integrals can be grouped into so-called families which are by definition disjoint. As we 
will see below this ordering is crucial in the construction of the coefficient functions. 

Let us start from the master integrals which have the most non-zero indices where only 
indices corresponding to the propagators but not to the numerators are counted. These 
integrals are maximal in the sense of the hierarchy defined above. For a maximal master 
integral F{nj), the corresponding coefficient function Cj(n) in Eq. (0) is replaced by the 
corresponding parametric integral (|7j) which results in a sum of terms as given in Eq. (fTTj) 
according to our agreement of choosing integration contours for the variables associated 
with the indices equal to one. We want to stress that only those integrations over Xj are 
non-trivial where the corresponding index of the master integral is zero. Experience shows 
that these remaining integration can be performed in terms of gamma functions, for the 
general integer values of the indices.^ At this point one applies Eq. (jH)) to normalize the 
resulting expression. 

Let us now suppose that we are dealing with a master integral which is not maximal, 
i.e., we have two master integrals which fulfill the hierarchy F{n2) < -^(lii)- This means 
that if = 1 we have also riu = 1. To construct an algorithm for the coefficient function 
C2(n) we start with the case of negative indices rij for those indices j where riu = 1 as in 
this case we have ci{n) = 0. Experience shows that the integrations for C2(n) result in 
ratios of gamma functions which in particular can be used to obtain 02(112) = 1- In a next 
step one considers the case nj > 0. Then the corresponding parametric representation (|7jl 
usually leads to integrals (cf. Eq. (fTT|) ) which cannot be evaluated in gamma functions.^ 
Thus at first sight it looks hopeless to achieve that the coefficient functions have to be 
rational functions of d (and eventually other dimensionful quantities). The way out of 
this problem is as follows: actually we have to look for an expression for the coefficient 
function C2(n) which is a linear combination of ci{n) and the representation ((Zj) for C2{n). 
Denoting the latter by 02(21) one has 



where the constant A is determined by the condition 02(211) = (^f- EQ- ®) and is thus 
given by 



At this stage it turns out that it is advantageous to use IBP applied to the parametric 
integrals of Eq. (fTTj) . Note that the complexity is significantly reduced as compared to 
the original integral simply because there are much less variables involved. The basic 
idea is to express any given parametric integral in terms of auxiliary (parametric) master 
integrals and expressions which are straightforwardly evaluated in gamma functions. The 
dependence on the new auxiliary master integrals has to drop out in order to guarantee 
that the coefficient functions are rational functions of d (and eventually the other dimen- 
sionful parameters). Such a cancellation turns to be a good check of the algorithm. This 

"^Typically, integrals involving the Euler beta function arise as we will see in the examples discussed 
in Section 131 

'"^However, sometimes this is possible as we will see in the examples below. 
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is similar to the cancellation of spurious poles when applying general prescriptions for 
asymptotic expansion of Feynman integrals in various limits of momenta and masses j3] . 
In our particular example of a hierarchy of two master integrals one presumably only 
needs one auxiliary parametric master integral. 

For a generic tower of master integrals the construction and cancellation of the aux- 
iliary master integrals goes along the same lines. When constructing the corresponding 
coefficient function one uses a linear combination of the coefficient under consideration 
and all coefficients of the higher master integrals. The constants are determined with the 
help of Eq. ^ . Explicit examples can be found in Section IHl 

Let us stress again that in general it is rather simple to compute the coefficient function 
of a maximal master integral as it has many indices equal to one. On the other hand, 
following the procedure described above, it is possible to either reduce the indices of the 
coefficient of a non-maximal master integral in order to make an evaluation in F functions 
possible, or relate it to the coefficient functions of higher master integrals. This is achied 
with the help of IBP applied to the parametric integrals over the x— parameters. 

At this point it may happen that the initial list of the master integrals obtained in 
Section ITT] is not sufficient and has to be extended by additional master integrals where 
some of the indices are negative. This necessity is connected to the type of solution of the 
recursive problem for the auxiliary parametric integrals and thus to the basic polynomial 

P^ix). 

We consider the construction of a coefficient functions as completed if for given indices 
an algorithm (usually realized on a computer) for its evaluation is provided. However, in 
simple situations, one can even obtain explicit expressions for the coefficients as functions 
of the indices. 

After we have obtained the list of the master integrals it is reasonable to check whether 
there are identical expressions among them, although they all have appeared in a different 
way. If some of the master integrals within one family are equal simply because the 
corresponding integrals are identical, e.g., due to some symmetry, it is natural to consider 
one master integral instead. The corresponding coefficient function is the sum of the 
coefficient functions of the original integrals. From the mathematical point of view, the 
best solution would be to achieve the absolute minimal number of the master integrals in 
Eq. (0). However, for practical purposes, this is not very important. 

Finally let us note that it is possible to identify an integral as master and then im- 
mediately construct an algorithm for the corresponding coefficient function, rather than 
strictly decompose the procedure into the two steps as described above. 
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(a) (b) (c) (d) 

Figure 1: Feynman diagrams for one- and two- loop massless propagator ((a) and (b)). 
At two-loop order there are two master integrals which are shown in (c) and (d). 

3 One- and two-loop examples 
3.1 One- loop propagator diagram 

Let us start with the simple example of a massless one-loop diagram with one external 
momentum (a so-called propagator diagram) which is given by (see also Fig. HJ^a)) 

where the usual prescriptions like fc^ = fc^ + iO, etc. are implied. According to [T], the 
transition to the corresponding vacuum problem reduces to adding a new propagator, 
l/(g^ — s), with some mass squared s. Furthermore one has h = 2. We want to consider 
the value of our diagram at some general point and are not interested in higher terms of 
the Taylor expansion in g^. Therefore we consider only the index 77,3 = 1 corresponding to 
the additional propagator. The integration contour for the corresponding variable X3 is 
taken as a Cauchy contour around the origin. Thus we write down the corresponding basic 
polynomial (jH)) and reduce it by setting 0:3 = 0. After that the resulting basic polynomial 
reads 

P(xi, X2) = {q^f - 2q\xi + X2) + (xi - . (15) 

The only possible candidate for a master integral is F(l, 1) because we obtain massless 
tadpoles when we put one of the indices to zero. We have 

,,.Fii,i). ,:.^/^(_,y/^-a r(2-^^/2)r^W2-i) ^^^^ 

The corresponding coefficient function is easily obtained from Eq. ((7j) with the help of the 
Cauchy theorem 



, (17) 

Xi=0 
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where A'^i is a (d-dependent) normalization factor introduced in order to fulfill Eq. ®. 
For our choice of P{xi,X2) it reads 

N,={ey"'. (18) 

It is straightforward to implement Eq. (|T7jl on a computer. The result agrees with 
the known one which in this case, of course, is trivially obtained by a straightforward 
Feynman parameterization. 



3.2 Two-loop propagator diagram 

The situation becomes slightly more complicated at two-loop level where the Feynman 
integral is given by (cf. Fig. Hl^b)) 

^(^) = / (A;2)ni(/2)n2p _ g)2]n3[(/ _ q)2^n,^^k - /)2]-5 " ^^^^ 

The transition to vacuum integrals is performed in analogy to the previous one-loop 
case by introducing one more propagator connected with the external invariant and 
h = 3 (see also the next example for more details). As before, we set the corresponding 
x-parameter to zero and obtain the following basic polynomial 

/ \ o 

P{Xi, . . . , X5) = —X1X2X3 + X2X3 + X2X3 -|- X1X4, — X1X2X4 — XiXsXi — X2XSX4 + XiX^ 

+ X1X2X5 - X2X3X5 - X1X4X5 + X3X4X5 + q'^[-XiX3 + X2X3 + XiXi 

- X2X4 + X1X5 + X2X5 + X3X5 + X4X5 - X5] + {q^Yx^ . (20) 

An analysis of the integrals with indices equal to zero leads to two master integrals shown 
in Fig. mc) and (d) 

/i = F(l,l,l,l,0) = (^vr'^/2)2(_^2)d-4^ Y2^d_2) ' 

/2 = F(0,1,1,0,1) = F(1,0,0,1,1) = -(^-'''n-<ir-'^-^^Y^^^ 

where the second one appears in two ways. As far as the coefficient functions are concerned 
the integration over those Xi with index "1" in the corresponding master integral is trivially 
obtained through differentiation of P(x). In the case of ci this leads to the following 
remaining integral 

Here and in the following a and (3 depend on the dimension d whereas k represent integer 
indices. On the other hand, for C2 one has to solve a two-dimensional integral over xi and 
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(a) (b) 

Figure 2: (a) Generic Feynman diagram needed for the one-loop heavy quark potentiaL 
There two master integrals: one corresponds to the very diagram with all indices equal 
to one and the second one is shown (b). 

X4 which is calculated as follows 



Jo 



13 



_ ^ r(ai + l)r(a4 + l)r(-ai -a^- (3-2) 

-WJ r(-/3) • ^^^^ 

Again a small computer program takes over the differentiation and integration parts and 
finally results to the two-loop integral 

F(ni, 77,2, na, n^, n^) = Ci(ni, 712, n^, n^, n^)Ii 

+ [C2(ni, ^2, ^4, ns) + 02(^2, rii, 714, n^, rig)] h , (24) 

in agreement with the usual approach based on IBP 1^. 

3.3 One-loop diagram for the heavy quark potential 

The integrals studied in this and the next subsection are useful in the context of non- 
relativistic QCD (NRQCD) where in addition to relativistic propagators (for gluons and 
light quarks) also non-relativistic ones (for heavy quarks) appear. The general one-loop 
integral is shown in Fig. |21 and has the form 

with v-k understood as v-k + iQ. Here q is the external momentum and f is a constant 
vector which fulfills the condition f -g = 0. 

When turning to the corresponding vacuum problem we consider, in addition to fc^, 
q-k and v-k, three extra kinematical invariants, which play the role of inverse 

auxiliary propagators. The latter are obtained after interpreting the expression in Eq. (j25|) 
as a three-point function with external momenta q and v and closing the three external 
lines in a new vertex. This formally leads to a three-loop vacuum problem and thus to 
h = ?>. 
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In the notation of Eq. ((T)) we have Di = k'^, D2 = {k — qY, = k-v + v'^, D4 = f^, D5 = 
and Dq = (g + f )^ which defines the matrix A according to Eq. The symmetrical 
extension of the corresponding inverse matrix leads to (after identifying the six invariants 
with xi, . . . , xq) 

Xi {Xi-X2 + X5)/2 X3-X4 \ 

{Xi-X2+X5)/2 X5 {-X4-X5 + Xq)/2\ (26) 

X3 - Xi {-Xi - X5 + Xq)/2 X4 J 

The basic polynomial P is obtained from the determinant of this matrix with the help 
of Eq. (jS}. Moreover, we shift the variables Xi by the corresponding effective masses: 
^3^X3 + v'^,Xi X4 + f^,X5 ^ X5 + q'^,XQ Xq + {q + f)^. Note, that as in the 
previous examples we are not interested in higher order Taylor coefficients of the additional 
kinematical invariants. Thus we can reduce the basic polynomial by setting X4 = x^ = 
Xq = and finally obtain 

P(xi, xs, X3) = {q^fv^ + (xi - X2f + 2g2 p ^ _ ^^f^ ^ (27) 
which leads to the two master integrals 

(_g2)d/2-5/2 0pr(5/2 _ dl2)T^{dl2 - 3/2) 



Ji = F(l,l,l) = -i-K^/''- 



V V{d-?,) 



u ^ 1.0) ^ .„-.(-,v/^-. r(2-^/2)r^W2-i) , 

Since for l\ all indices are one the coefficient function c\ is obtained in analogy to 
the one of Eq. (fT7j) . In the case of C2 the indices n\ and n2 are treated with the help of 
differentiations which leads to the following type of integral for X3 



r(3) 



{k3, a) = J" dx3 (a^ - ^l)" . (29) 



Here /cs is an integer but a depends on d. {k, a) can be interpreted as a principal value 
integral leading to 



/f(fc,a) = <|ry r(« + fc/2 + 3/2) ^--^-v__ ^^^^ 

for k odd 
The final result 

F(ni, n2, ris) = Ci(ni, ^2, ns)/! + C2(ni, ^2, n3)/2 , (31) 
agrees with the explicit analytical result (see, e.g., Ref. flSj)- 




a+fc/2+1/2 r(A;/2 + l/2)r(a + 1) , , 

for k even 
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(a) (b) 
Figure 3: Feynman diagrams corresponding to case A (a) and case B (b). 

3.4 Two-loop diagrams for the heavy quark potential 

At two-loop order two classes of Feynman integrals are needed which we will refer to as 
case A and case B 

~ J J (A;2)"l(/2)"2[(A; - g)2]n3[(/ _ g)2]n4p _ /) 2] "5 . A;)n6 . /)nr ' ^^2) 
Fein) - J J (p)„,(p)n2p _ g)2]n3[(/ _ g)2]n4p _ /) 2] "5 . A;)"6 . ( A; - /)]"r ' ^^^^ 

They are shown in Fig. El The corresponding basic polynomials read 

Pa(xi, . . . , X7) = -[X2X6 - X4XQ + (-Xi + X-i)x7]'^ + f ^{XiX4 + X3{xl + X2{x3 - X4 - X5) 
+ X4X5) - Xi[x2{Xs + X4- X5) + X4(X3 - X4 + X5)]} 

+ (g^)^[t;^X5 - {xq - Xjf] + q'^{v'^[{x3 + X4- X5)X5 + X2{X3 -X4 + X5) 
+ Xi(-X3 + X4 + X5)] + 2[X2X6(-X6 + X7) + X4X6{-Xq + Xj) 

+ Xt{xiXq + X3X6 — 2x5X5 — X1X7 — X3X7)]} , (34) 

Pb(xi, . . . , X7) = Pa{xi, X2, X3, X4, X5, X6, X6 - X7) . (35) 

The application of the procedures described in Section|21to case A leads to the following 
families of master integrals which are also pictured in Fig. |3 As far as the notation is 
concerned the first index labels the different master integrals. In case the master integrals 
are equal we introduce a second index for further specification. An additional subscript 
"a" is added for those master integrals where one index is "—1". 

• Family Al. One has four master integrals which obey the following hierarchy: 
h > {/2i,/22} > ^3 with 

/i = Fa(1,1,1,1,0,1,1), 
/2i = F^(l,l,l,l,0,0,l), 

J22 = F^(1,1,1,1,0,1,0), 

/3 = F^(1,1,1,1,0,0,0). (36) 

• Family A2. One has four master integrals which obey the following hierarchy: 

hi > {^7i,^8i} > hi with 



'51 



F.fl, 0,0,1,1, 1,1) 
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Family Al: 



Family A2: 




-^51 hi hi hi 



Family A4: 




hi 

Figure 4: Feynman diagrams corresponding to the master integrals of case A. In addition 
to Iqi there is also a master integral (/eia) containing an irreducible numerator (see text). 

Jn=F^(l,0,0,l,l,0,l), 
J8i=Fa(1,0,0,1,1,1,0), 

/4i=F^(l,0,0,l,l,0,0). (37) 

• Family A3 is symmetrical to Family A2 with respect to 1 ^ 2, 3 ^ 4, 6 ^ 7. It 
contains the master integrals 1 52, h2, h2 and 142. 

• Family A4 contains the master integrals 

/6i = Fa(0,1,0,1,1,1,0), 

/61a = i^A(0, 1,0, 1,1, 1,-1). (38) 

• Family A5 is symmetrical to Family A4 with respect to 1 2, 3 ^ 4, 6 7. It 
contains the master integrals Iq2 and /62a- 

As expected, the master integrals of the massless scalar two-loop diagram discussed in 
Section also appear in the above list. 

As has already become clear from the examples discussed so far, one expects the 
appearance of complicated expressions for the coefficient functions of the simple master 
integrals. Indeed, in the case of the coefficient function ci six out of seven indices can 
be treated with the help of differentiations and the remaining one-dimensional integral 
is easily solved with the help of ij^^ given in Eq. ()22|) . The situation is similar for C22 
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(and C21 which can be obtained by exploiting the symmetry) where the remaining two-fold 
integration over X7 and X5 is performed with the help of 



J a 



^ ' {k-r)\r\ r(ai + a2 + r + 2) ^ ' 

and Eqs. and (IHUj). respectively. In principle, instead of performing the x^ integration 
one could also introduce an auxiliary master integral as a result of a two-dimensional 
(xs and the dimension d) recursion problem. This master integral would cancel after 
considering the proper linear combination with ci as described around Eq. p2|l . 

For the coefficient C3 there are three non-trivial integrations (0:5 , xg, xj) left. In case 
one of the indices 725, rig or rtj is less or equal to zero one can use various combinations 
of the auxiliary integrals I^^ {i = 1, ... ,4) listed above. Thereby it is advantageous to 
perform the integration corresponding to the negative index first. If, on the contrary, ns, 

and nj are positive an immediate integration seems not to be possible. However, from 
the corresponding three-parametric integral representation it is simple to derive recurrence 
relations which shift at least one of the indices to zero, eventually at the cost of increasing 
the dimension. The latter does not constitute a problem since the whole formulation of 
our procedure is in d dimensions. Thus, also in this case the integration can be performed 
in terms of F functions. In principle one could be forced to introduce three auxiliary 
master integrals and build the proper linear combinations with ci, C21 and C22. However, 
it turns out that the corresponding constants, 03(27,^) {i = 1,2,3) are zero. 

For the coefficient function C51 only two non-trivial integrations over X2 and 0:3 are 
involved which can be performed with the help of Eq. ()23|) . 

The situation is similar for C71. It is convenient to perform in a first step the Xq 
integration with the help of ij^^ of Eq. (j29p . For the remaining integration over X2 and 
X3 the formula 

II {k, a2, as, (3) = dx2 dx^ x'^^x'^^ [cf' + X2 + X3 j [(f + X2 j 

i2+a3+/3+fe+2 F(a2 + l)r(a3 + l)F(-a2 - a^- [3 -k-2) 



x /<-°'/:'' . (40) 

r(— 0:3 — p — k — 1) 

turns out to be very useful. In principle, the complete construction of C71 requires the 
proper linear combination with C51. However, since we have 071(215^) = this is not neces- 
sary. For cgi no separate calculation is necessary as the symmetry of the basic polynomial 
can be exploited and one has C8i(ni, n2, n3, n4, ns, ng, Ti?) = 071(114, n3, n2, ni, 7x5, 117, ng). 

The most complicated coefficient function is certainly C41 since there are four non- 
trivial integrations over X2, X3, xq and X7 left. If ng or rii are less than or equal to zero 
the integrations can be performed in terms of F functions with the help of the formulae 
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provided above. However, for > 1 and > 1 this is not possible. In this case the idea 
is to use IBP in order to reduce the four-parametric integral representation 

7-A,aux/ N f f \r> / \^z-ni dx2 da;3 da^e dxy 

J J X'2 Xq Xij 

(with z = {d — h — l)/2 = d/2 — 5/2) to the auxiliary master integral 1, 1, 1, 0). 

P41 is obtained from P4 by setting xi, X4 and X5 to zero. We should stress that the 
corresponding recurrence procedure is significantly simpler than the original one which 
involves seven denominators. Furthermore, if during the recursion either tiq or n-i be- 
comes negative the corresponding expressions can immediately be expressed in terms of F 
functions. The five IBP relations which are useful for the reduction to /^'^"^^(l, 1,1,1,0) 
can be obtained by either differentiating the integrand with respect to Xi {i = 2, 3, 6, 7) or 
by writing P/f"'^ = Pi\^'^~^Pii and inserting the exphcit result for the last factor. The 
proper combination of these relations leads to new ones which allows the following steps 
to be performed in an automatic way: 

1. Reduce uq and n-j to one. 

2. Reduce n2, na > to n2, < 0. 

3. Use IBP recurrence relations to obtain 712 = ris. 

4. Reduce n2 — < to n2 — ris — 0. 

5. Adjust the dimension, i.e., reduce rid to zero. 

A simple relation transforms /^''^"''(O, 0, 1, 1, 0) to /^''^"^(l, 1, 1, 1, 0). At this point one 
constructs the final coefficient function C41 by considering the linear combination with C51, 
C71 and Cgi. Since one has 041(2171) = C4i(n8i) = we are left with 



where 



ifei) 



C4i{n) = C4i(zi) - C4i(2l5i)c5i(zi) , (42) 



1 4(rf - 3) {3d - 14) {3d - 10) {3d - 8) 



4i\mij- 02 



(rf-4)2(3(i- 13)(3(i- 11) 

^3<i -'13^(3! 11) (^^) 

In this combination the auxiliary master integral /^'^"^(l, 1, 1, 1, 0) cancels and 041(77,) is 
a rational function in d. 

The master integral /ei forms a family by its own. However, as the polynomial Pei 
is quadratic in X7 and thus the corresponding recurrence relation shifts nr only in steps 
of two it is necessary to introduce in addition the master integral /gia where 777 = — 1. 
The practical implementation within our method does not need this recurrence relation. 
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Family Bl: 



Family B2: 



Family B3: 



Family B4: 




Figure 5: Feynman diagrams corresponding to the master integrals of case B. In addi- 
tion to Iqi {i = 3, . . . , 7) there are also master integrals (/eia) containing an irreducible 
numerators (see text). 

however, feels the necessity of introducing /gia. The very calculation of the coefficient 
function is identical for Iqi and Jgia- For < it can be done in terms of F functions 
with the integration order X3, xi, xj. On the other hand, for 77,3 > a simple one-step 
relation reduces to zero. 

Let us now come to case B. As one can see from Eq. (jH^j) the basic polynomial is quite 
similar to the one of case A which can be used while computing the coefficient functions. 
However, the symmetry can only be exploited if 727 < as for riy > the factor [x^ — Xj) 
would appear in the denominator. 

Altogether there are four families which, however, show a more complicated structure 
than in case A. More precisely one has 
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Family Bl. There are twelve master integrals which obey the hierarchies If > 
{/f , I22} > h and /f > /f > {hi, /em} {i = 3, 4, 5, 6) and are given by 



Fb(1,1,1,1,0,1,1 
.2 -F5(l,l,l,l,0,0,l 
722 = Fb(1,1,1,1,0,1,0 
73 = Fb(1,1,1,1,0,0,0 

763 = i^s(l,l,l,0,0,0,l 

764 = Fb(1,1,0,1,0,0,1 

765 = Fb(1,0,1,1,0,0,1 

766 = Fb(0,1,1,1,0,0,1 

Furthermore there are master integrals with uq = —1 



/63a = 1,1,0, 0,-1,1), 

/64a = i^B(l,l, 0,1, 0,-1,1), 
/65a = i^B(l,0, 1,1, 0,-1,1), 
/66a = /^b(0,1, 1,1, 0,-1,1). 



(44) 



(45) 



Family B2. One has four master integrals which obey the following hierarchy: /g > 
{/82,/8i} > hi with 

/9 = Fb(1,0,0,1,1,1,1), 
/82 = 7^^(1,0,0,1,1,0,1), 
781 =FB(l,0,0,hLLO), 

/4i = Fb(1,0,0,1,1,0,0). (46) 

Family B3. Similarly to Family B2 one has four master integrals obeying the hier- 
archy hs > {hs, h2} > h2 with 

/53 = 7^b(0,1,L0,1,1,1), 
783 = 7^b(0,LL0.L0,1), 
/72 = Fs(0,l,l,0,l,l,0), 

/42 = Fs(0,l,l,0,l,0,0). (47) 



Family B4 consists of the two master integrals 

/67 = 7^b(0,1,0,LL1,0), 

/67a = /^b(0, 1,0, 1,1, 1,-1), 



(48) 



which is similar to the Famihes A4 and A5 of case A. 
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The construction of the coefficient functions cf , cf and C22 from the family Bl pro- 
ceeds along the same lines as in case A. In the case of C3 there is a slight complication 
as, in contrast to case A, C3{ni) 7^ 0. As a consequence an auxiliary master integral^, 
jB,aux^g^ 1, 1, 0) has to be introduced which is only canceled after considering the proper 
linear combination with ci. The reduction to I^'^^^{0, 1, 1, 0) is straightforward. 

Family Bl has four more members, Jes, Iqa, /es and fee, which belong to the four 
hierarchies > I2 > Im {i = 3, 4, 5, 6). Thus, in order to obtain the coefficient functions 
cgj one has to consider the linear combination 

ce^ = 4 - cUnf)cf{n) - 4(nf )cf (n) . (49) 

Let us in the following restrict the discussion to Cgs since the results for the other three 
coefficients can be obtained by exploiting the symmetry. Cgg is given by an integral 
representation of the form 

f f \D / \^z-nd dX4dX5dX6 , - 

with 

^63 = {q^fv'^x^ + q^v"^ (xiX^ - xfj - Aq^x^xl - x\xl . (51) 

For 77-4 < 0, where we have cf (ra) = 0^(22) = 0, the integrals in Eq. ()50|) can be solved 
analytically in the order 0:4, X5, Xq using Eq. for X4, the formula 

I^^\a,p) = J^dxx" {x + af 

L+i r(i + ^)r(-a-/3-i) 

n-p) ' ^^^^ 

for and Eq. extended to non- integer for xq. For 77,4 > two auxiliary master 
integrals, /^''^"''(l, 0, 0, 0) and /^'^"""(l, 0, 1, 0), have to be introduced where the reduction 
of Eq. proceeds s follows 

1. Reduce to one. 

2. Reduce to zero. 

3. The reduction of uq can only be performed in steps of two. Thus one end up with 
riQ = or riQ = —1. 

4. Adjust the dimension, i.e., reduce to zero. 



^The definition of /f '^"'^(77.5, ng, ny, Ud) is in analogy to Eq. (141(1 . 
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The corresponding recurrence relations are easily derived from Eq. ()51|). It is interesting 
to note that in Eq. ()49|) the master integral /f^'^"^(l, 0, 1, 0) is canceled from cf and 
0, 0, 0) from cf . Please note that due to the structure of Eq. (j5T| in addition to 
Jes also a master integral with uq = —1, /gsa, has to be introduced which, however, has 
the same coefficient function as /gs. We want to mention that for cga and cgs the master 
integrals Iq and Jga are needed while for cq^ and Cqq the integrals Iq and If? are necessary. 

As far as the families B2, B3 and B4 are concerned the strategy to construct the 
coefficient function goes along the same lines as for the families A2, A3 and A4. 

For completeness let us in the following list all occurring master integrals. They 
have been obtained with the help of the program package developed for the calculation 
performed in Ref. JHI where conventional recurrence relations have been applied. In the 

results the overall factor (j.tt'^^'^^ Q~^^ is omitted and Q = a/— ^ has been introduced. 
Moreover the standard factor e"^''''^^ is pulled out in the formulae where an expansion in 
e is performed. We obtain 



h 
h 
h 



71 T\5/2-d/2)T\d/2-3/2) 

Qh^ r2(d-3) ' 

^r{2 - d/2)T{h/2 - d/2)T'^{d/2 - l)T'^{d/2 - 3/2) 



Qv r(rf-3)r(rf-2) 

V\2~d/2)T\d/2-l) 



T^{d-2) 

h = -Q 



2r(3-rf)r3(c//2- 1) 



r(3rf/2 - 3) 



h 
h 



TT 
~i2 



2^2 , 7 



-^-4 + (-24+-7r2)e + 0(e2) 



v^g 2^-^T{3 - d)T{7/2 - d)T{d/2 - l)T\d - 5/2) 



ha = -\^Q 



r(2-d/2)r(2rf-5) 

2^-2^2(3 _ d)T{d/2 - l)T\d - 2) 



T{3/2-d/2)T{2d-A) 
nQ r(7/2 - d)T\d/2 - l)T{d/2 - 3/2)T{d - 5/2) 



V r(rf-2)r(3rf/2-4) 

h = l7, 

h = h, 

1 



1 2 ^' 



-'2 



-41n2 + e (^^TT^ - 161n2-41n22) + C(e^) 



Qv 

lia = -ha ■ (53) 

Observe that h = h and h = h although the corresponding integrands are certainly 
different. These equalities can be immediately be seen by a simple change of variables. 
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Because of I7 = Is we have in both cases one master integral less and end up with 
eight master integrals, Ji, . . . , J7 and Jga, in case A while in case B ten master integrals 
contribute, I2, . . . , Ij, Ig, Igay ^^nd Please note that only two of them are not known 
in terms of F functions. In the above list their result is given as an expansion in e up to 
terms of order e^. 

A realistic diagram (in particular for arbitrary gauge parameter ^) does not imme- 
diately lead to the integrals Fa and Fb defined in Eqs. (j22I) and (jHHI) . However, it is 
straightforward to map the expressions to the function (see also the appendix of Ref. [T7j ) 

^-'^ J J (p)ni(/2)n2p _ g)2]"3[(/ _ g)2]n4[(^ _ lfYB[vkY^{vlY^[v{k - ' 

(54) 

with the obvious identities 

J(ni,n2,n3,n4,n5,n6,n7,0) = F^di), J(rai, ng, ng, 714, ng, ng, 0, ng) = FB{n). 

The reduction to Fa and Fb is achieved with the help of the relation 

J{lL) = J(ni,n2,n3,n4,n5,n6 + l,nj- l,ns) + J(ni, na, ng, n4, ng, ne + I,n7,n8 - 1) , 

(55) 

where either rii or ng are reduced to zero from positive values. In the case of negative 
or ng such a reduction is similar. 

The method described above has been implemented in a Mathematica package. The 
typical runtime of an integral contributing to the static potential or to the l/lmgr"^) 
corrections amounts to a few seconds. This extends to the order of a minute for integrals 
which occur in the calculation for a general gauge parameter where the exponents of the 
propagators can be significantly higher. 

Let us at the end of this section present some explicit results which illustrate the 
functionality of our method. For this reason we list next to the results for FA{n) and 
Fb{r) also the (non-zero) expressions for the individual coefficient functions. 

For n = (2, 2, 1, 1, 1, 1, 1) we get^ 

1 / 2 4 , 16 368 , s s\ 

F^(2, 2, 1, 1, 1, 1, 1) = + - y + - 8^(3) + 0(6)) , (56) 

with the coefficients 

2(d - 5)(d - 4) 8(d - 5)(d - 3f 

ci = 7: , C3 - 



g6 ' = (rf-4)gV ' 

_ -3{d - 3){3d - 16)(3rf - U){3d - 10) (3c/ - 8){5d^ - 93^^ + 588rf - 1264) 
~ ~ [d -9){d- 8){d -7){d- 6)2(rf - 4) V°^2 

-3{3d - 17)(3rf - 13)(3rf - 11) 



C51 — C52 



{d~9)id-7)q^ 



_ -32{2d - 13){2d - ll){2d - 9){2d - 7){2d - 5) 
~ ~ id-9){d- 7)i d -6)id- 4")^ ' ^ ^ 



^The same standard overall factors are pulled out as in Eqs. H53|l . 
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where cei and cq2 get multiplied by I^a- In the case B the results read 

Fsi2, 2, 1, 1, 1, 1,1) -r^. (-^ + + ^ + ^^(3) + ^(^)) ' (^S) 



. _ 2(d-5)(d-4) _ -4(d -5){d- 3Y 

Cl — a , C3 



C41 = 
C42 = 
C53 = 

C63 — C64 = 
C65 = C66 = 
C67 



g6 ' " (d-4)gV ' 

3((i - 3)(3d - 16)(3c^ - 14)(3ci - 10)(3ci - 8){7d^ - lllcP + 654(i - 1232) 

{d -9){d- 8){d -7){d- 6)2(rf - 4)V°t'^ 
-6((i - 3) (3(i - 16) (3rf - 14) {3d - 10) {3d - 8) (ci^ - Ud"^ + 33d + 16) 

{d -9){d- 8){d -7){d- 6)2(d - 4)2gio^;2 ' 
-3(3d - 17)(3(i - 13)(3d - 11) 



{d-9)(d-7)q^ 

4(2d - 7) (2d - 5)(15d^ - 304^3 + 2240^^ - 7093^ + 8118) 
{d -9){d- 7){d -6)(d- 4)gi%2 

4(2(/ - 7) {2d - r)){d- - 17 d + 55) 
{d - 7){d - 4)q^%^ ' 
-32(2d - 13)(2d - ll)(2d - 9){2d - 7)(2d - 5) 



{d-9){d-7){d-6){d-4)q^0v 



2 



-3(3d- 17)(3ci- 13)(3d- 11) 
^ {d-9){d-7)q^ ' ^'^^ 

where cgs and ces get multiplied by Iqu and C64, Cee and cer by As a second example 
let us consider n — (1, 1, 2, 1, 1, —1, 1) which gives 

Fa{1, 1, 2, 1, 1, -1, 1) = ^ (-1 + ^ - 2C(3) + 0(e)) , (60) 
2{d-3) -3{3d-10){3d-8){d'^ -5d + 2) 



C42 



(d-4)g4' 2(d-6)(d-5)((i-4)V 
3(d - 5)(ci - 2)(3d - 10)(3d - 8) 
2(d-6)(d-4)V ' 



4(2d-9)(2ci-7)(2d-5) 
= {d-5){d-A)q^ ' ^^'^ 

where cq2 gets multiplied by /ea- 

Fb(1, 1, 2, 1, 1, -1, 1) = ^ (-^ + ^ + ^(e)) ' (62) 

_ (d-5){d-3) _ -3(3d - 10) {3d - 8) (d^ -9d+ 22) 

~ (d - 6)g4 ' - 2(d-6)2(d-5)(d-4)g6 ' 

22 



_ 3{3d - 10)(3rf - 8)(d2 - lid + 26) _ {2d - ll)(2rf - 7){2d - 5) 

~ 2(rf-6)2(c/-4)g6 ' ~ (rf-6)(rf-5)g6 

(2ci-7)(2rf-5) (2rf-7)2(2d-5) 



•^64 — — n — T^vn — r\ fi ' "^65 
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(d-6)(rf-5)g6 ' °^ ((i_6)(rf-5)g' 
_ -(2rf-7)(2rf-5)(4rf-19) 

- {d-6){d-5)q^ ' ^^^^ 

where again cea and ces get multiplied by I^a and C64 and cee by 



4 Conclusion 



In this paper we described an algorithmic procedure to identify the master integrals for 
a given class of Feynman diagram. Furthermore details are given for the computation 
of the corresponding coefficient functions which are formally defined in Eq. ((7j). The 
general procedure was applied to useful one- and two-loop integrals which often appear 
in practical applications. In particular we discussed in great detail the integrals needed 
for the calculation of the two-loop heavy quark potential. We do hope that, with our 
instructions, the reader will be able to apply Eq. ((Tj) to various recursion problems for 
Feynman integrals. 

Let us finally point out that there is also another approach to use Eq. d?)) which is 
based on an expansion of some parametric representations (originating from Eq. ((7j)) in the 
limit of large dimension d ^21- There the crucial assumption is the rational dependence 
of the coefficients Cj on d. No description of this procedure has been published up to now. 
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